Photonic Crystal Nanobeam Cavity Strongly Coupled to the Feeding Waveguide 
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A deterministic design of an ultrahigh Q, wavelength scale mode volume photonic crystal 
nanobeam cavity is proposed and experimentally demonstrated. Using this approach, cavities with 
Q> 10® and on-resonance transmission T>90% are designed. The devices fabricated in Si and 
capped with low-index polymer, have Q=80,000 and T=73%. This is, to the best of our knowledge, 
the highest transmission measured in deterministically designed, wavelength scale high Q cavities. 
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Photonic crystal (PhC)[l|[l| cavities, with quality(Q) 
factors over million and wavelength scale mode vol- 
umes, are widely applied in the fields that range from 
quantum information processing and nonlinear optics, 
to biomedical sensing. The fomier space analysis Q [1] , 
the multi-pole cancelation^, and the mode matching 
mechanism'^ [7] [s"] have been developed to explain the 
origin of high Qs. The design of PhC cavities, how- 
ever, is typically based on extensive parameter search 
and optimization [o'l-fisl], also known as intuitive design. 
The large computational cost, in particular the compu- 
tation time, needed to perform the simulation of high- 
Q cavities make this trial based method inefficient. In- 
verse engineering design, in which the physical structure 
is optimized by constructing specific target functions and 
constraints, was also proposed 14 1 [l^ . A design recipe 
based on the desired field distribution is proposed in [1^ . 
In this letter, we propose and experimentally demon- 
strate a deterministic method to design an ultrahigh 
Q, sub- wavelength scale mode volume, PhC nanobeam 
cavity(FigurelT]) that is strongly coupled to the feeding 
waveguide(i.e. near unity on resonance transmission). 
The design approach is deterministic in the sense that it 
does not involve any trial-based hole shifting, re-sizing 
and overall cavity re-scaling to ensure ultra-high Q cav- 
ity. Moreover, the final cavity resonance has less than 
2% deviation from a predetermined frequency. Our de- 
sign method requires only computationally inexpensive, 
photonic band calculations (e.g. using plane wave expan- 
sion method), and is simple to implement. 

The Q factor of a PhC nanobeam cavity can be max- 
imized by reducing the out-of plane scattering((3sc) due 
to the coupling to the radiation modes. As shown 
previously [3] [16], scattered power (Psc) can be expressed 
as an integral of spatial fourier frequencies within a 
light cone, calculated over the surface above the cavity: 
P^c ^ U^^tconedH\FT{H.)\' + \FT{E^)?). The in- 
tegral is minimized when major fourier components are 
tightly localized (in k-space) at the edge of the first Bril- 
lioun zonejl]. We start by considering the ideal field dis- 
tribution on this surface which would minimize -Psc ■ 
general property of these nanobeam cavities is that it 
consists of the waveguide region of length L, that sup- 
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FIG. 1: (a)Schematic of nanobeam PhC with conventional 
Bragg mirror. (b)SEM image of a silicon nanobeam PhC with 
"modulated Bragg mirror" . (c)&(d) logjg|FT(_ff)| for a model 
cavity with different cavity length L (normalized to the period) 
and two different types of mirror that support: c)exponential 
attenuation with k — 0.6 and d)Gaussian attenuation with 
a — 0.1. (e)The portion of the fourier components that are 
inside the lightline assuming an index propotion of 2.5 be- 
tween the waveguide and the surroundings. 



ports propagating modes, surrounded by infinitely long 
Bragg mirror on each side(FigurelT^). Without the loss 
of generality, we consider the TE-like cavity mode with 
Hz as a major field component. In the case of conven- 
tional periodic Bragg mirror, evanescent field inside the 
mirror can be expressed as sin(/3Bragga;) exp(— kx), where 
K is attenuation constant. The cavity field inside the 
waveguide region can be represented as sm{j3„gx). As 
mentioned above, scattering loss decreases in mirror sec- 
tion when /Ssragg = 7r/a, while phase matching between 
mirror and waveguide 0, /^Bragg = /3wg, minimizes the 
scattering loss at cavity-mirror interface. The spatial 
fourier transform of such cavity field is approximately 
a Lorentzian in the vicinity of n/a. As proposed in 
spatial frequencies within light-cone can be minimized 
(Q maximized) if the field attenuation inside the mirror 
has a Gaussian shape sin(/3Bragga;) exp(— cra;^). In Fig- 
ure [T](c)&(d), the fourier space spectrums of the cavity 
modes with exponential and with Gaussian attenuation 
are shown respectively. The fraction of the energy associ- 
ated with spatial harmonics within the light-cone, rj, for 
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both cases is shown in Figure[Tfe). It can be seen that 
cavity with the Gaussian attenuation has more than an 
order of magnitude smaller ry. 

The preferred Gaussian attenuation can be obtained 
by making k a linearly increasing function of the posi- 
tion within the mirror(K — ax). We name such mirror 
"modulated Bragg mirror" . At the same time, to prevent 
scattering, the oscillating part of the field, sin(/3inoda;), 
should have constant /3niod throughout the mirror. This 
condition is satisfied only if each segment of modulated 
Bragg mirror has the same length ("periodicity") a, and 
if the operating frequency is kept inside the bandgap of 
each segment of the modulated Bragg mirror. In the case 
of dielectric-mode cavity, this can be achieved by lower- 
ing the band-edges of each modulated mirror segment in 
a way that linearly increases the mirror strength. More- 
over, the optimal cavity length for this type of cavity is 
L = from Figurelljb), which simultaneously minimizes 
the mode volume (V) of the cavity and thus increases 
Q/V. It is worth noting that L/a — 3 in Figuredfb) 
corresponds to the so-called L3 cavity :9J. 

One way to make " modulated Bragg mirror" is to lin- 
early decrease the filling factor(FF) in nanobeam cavity, 
which, as we show later, correspond to a linear increase in 
mirror strength around the cavity center. FF is defined 
as the ratio between the hole area and the area of the 
unit cell. Other types of modulation may exist that en- 
able a linearly increasing mirror strength. In such struc- 
ture, cavity resonance frequency is expected to be very 
close but slightly smaller than the dielectric band-edge of 
the the central segment. This is due to the reduction of 
the hole size in the modulated Bragg mirror, and can be 
estimated using perturbation theory. The difference de- 
creases as the number of modulated mirror segments in- 
creases. To summarize the design principles we obtained 
so far: (i) zero cavity length {L = 0), (ii) the length of 
each segment be the same (period=a) result in constant 
phase velocity at ir/a and (iii) a modulated Bragg mirror 
results in Gaussian shaped field attenuation. 

We demonstrate the power of the recipe by designing 
an ultra-high Q and small mode volume PhC nanobeam 
cavity that operates at 1.525/ir7i(196.6THz) in a realis- 
tic geometry. We assume that the nanobeam is made 
with silicon-on-insulator material with 220nm thick Si 
device layer (constrained by our SOI wafer properties), 
and capped with silica: nanobeam is made of silicon 
(n=3.46), while holes, bottom and top cladding are made 
of material with (n=1.45) (air-holes backfilled with sil- 
ica). Our design approach is as follows: (i) We choose the 
period a by selecting ricff to be between ngi = 3.46 and 
JT-ciad = 1.45. rieff = 2.5 is a good compromise, resulting 
in a = Ao/2noff ~ 300nm. We note, however, that any 
UeS (and thus any periodicity) which opens a band gap 
can be used as a starting point to realize a high Q cavity, 
assuming that sufficiently slow modulation is used, (ii) 
Next, we choose the width of the nanobeam to be as wide 
as possible in order to push the mode away from the light- 
line, while still being single mode. Band diagram simu- 



lation shows that 700nm is an optimal choice, that keeps 
the second order mode (of the same symmetry as the fun- 
damental one) at the edge of the banggap. (iii) Next, we 
find the proper FF which produces a dielectric band edge 
at 1.525/im. With Uctf — 2.5, FF can be estimated us- 
FF)/3.462- 



17[ , which gives 



ing = (l-FF)/3.462-|-FF/1.452 

0.19. Then numerical band diagram simulation shows 
the actual FF = 0.15. Figurel2Kb)&(c) shows the band- 
diagram and mirror strength at several FFs. The linear 
region between FF=0.15 and FF=0.09 can be used to 
construct the "modulated Bragg mirror", (iv) Finally, 
the nanobeam cavity needs to be strongly coupled to the 
feeding waveguide in order to achieve a large transmis- 
sion efficiency when probing the cavity. The transmission 
at the cavity resonance can be written as Q^g/OtotaiE^- 
Therefore, the cavity Q should be limited by Qwg- The 
linear decrease of the FF from 0.15 to provides a natu- 
ral way to achieve this. More sophisticated couplers will 
be discussed elsewhere. In the proposed design method, 
(i) - (iv) , all cavity parameters are determined using fast 
band diagram calculations, only. We note, that the most 
critical part of our method is that the modification of the 
periodic photonic crystal is achieved by keeping the peri- 
odicity a constant, modulating the filing fraction and us- 
ing zero cavity length. This approach preserves the phase 
velocity of each segment and is essential for realization 
of high-Q cavities. Our design strategy has an additional 
important advantage over other types of photonic crystal 
nanobeam cavitiesfa- 12|,[l^-[21|: it provides a natural 



way to efficiently over-couple the cavity to the feeding 
waveguide. Again, we stress that these high Qs are ob- 
tained by-design, performing only simple band-diagram 
calculations, and no additional parameter-search trial- 
based method was performed. This significantly reduces 
the computational time to a few minutes, while a full 
FDTD simulation takes more than 24 hours using a grid 
with 64 processors. 



TABLE I: 3D simulation result of waveguide-coupled cavity 
with 40 and 50 modulated grating sections. The FF is 
changed from 0.15 at the center of the cavity to 0. Qwg and 
Qsc refer to the coupling Q factor to the feeding waveguide 
and the scattering Q to the radiation loss. Vcs is the mode 
volume normalized by (Ao/3.46)^. 



type 


X{fim) 


Qsc 


Qwg 


Qtotal 


Kff 


Trans. 


mod40 


1.552 


2.2E7 


1.3E6 


1.2E6 


1.1 


0.91 


mod50 


1.549 


2.5E7 


2.2E7 


1.2E7 


1.2 


0.53 



To verify our designs, 3D FDTD modeling is used to 
study the cavity with above mentioned parameters. The 
results are summarized in TableU for two cavities with 
40 and 50 modulated mirror segments on each side. It 
can be seen that cavities feature ultra-high Qs (both Qsc 
and Qwg)- At the same time, high transmission, partic- 
ularly in cavity "mod40", is obtained. In principle, Qs 
can be arbitrarily high by applying sufficient slow mod- 
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FIG. 2: (a)Schematic of modulated nanobeam cavity. (b)TE 
banddiagram of the above cavity with FF=0.15 and FF=0.07. 
The resonance of cavity "mod40" is about 1.5% lower than the 
dielectric bandedge of the central section with FF=0.15, due 
to modulation. (c)Mirror strength at different FFs. A linear 
increase can be seen at cavity center. (d)Simulated Ey pro- 
file at the middle of the naobeam cavity for the fundamental 
mode(Q=2.7E6) and second mode(Q=8E4). (e)Experiment 
transmission spectrum of the above mod40. Inset is the sim- 
ulation result of the bandedge modes of the current structure. 
Due to very large photon life time of our ultra-high Q cavity 
(^photon = Q/oj ~ Ins), it becomes nearly impossible to model 
transmission through the cavity using 3D FDTD. Hence the 
high Q cavity mode does not appear in the simulated spec- 
trum. (f)Zoom in of the transmitted signal of the fundamental 
mode at different input power levels. Dots are experimental 
data and lines are fitted curve using nonlinear equation 
The power levels in the legend is the power at the output of 
the fiber tip. Power coupled to the cavity is smaller due to 
the spot-size converter. 



ulation, at the cost of larger mode volumes. However, in 
practical case, cavity Q is limited by material losses and 
fabrication precisions, and therefore, Q of tens of millions 
is sufficient. 

To experimentally verify our designs, we fabricated 
waveguide-coupled cavity "mod40" using a silicon-on- 
insulator(SOI) wafer with a 220nm device layer on a 2um 
buried oxide, using the fabrication procedure described in 
our previous work^lSj. Spot size converter (23i. that con- 
sists of 2/xm X 2/im(cross section) polymer pad(npad = 
1.58) was used to couple light in and out of the cavities. 
The devices were embedded in polymer with refractive 
index ndad — 1.34. The reduction of refractive index of 
top cladding from the simulated one {n = 1.45) slightly 



affects our cavity design. We modeled devices with poly- 
mer top cladding and found that the cavity resonance 
is shifted to shorter wavelength (Amodci — 1517nm), the 
cavity Qtotai is 2.72E6, and the the on- resonance trans- 
mission is r = 75%. 

We characterized the device using a tunable laser 
source(1470-1580nm) and a tapered fiber tip(2.5/im spot 
diameter) to couple light in and out of the polymer 
waveguide. A polarizer is placed at the output to fil- 
ter out the TM-like mode. Figurej^l^e) is the experiment 
transmission spectrum through one of the resonators. 
The signal is normalized by the band-edge modes, which 
has a unity transmission as verified by 3D FDTD sim- 
ulations (inset). A nonlinear bistability is observed and 
Figure|2{f) further shows the zoom in of the fundamen- 
tal mode at different input power levels. We fitted the 
experimental data using the following expression typical 
for the nonlinear bistability 
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wg 



1 + {Pout/Po - 2(A - Ao)/7o)2 



(1) 



Po = 3K(3totalQwg(w/(2nc)) 

2^(3) is 

the characterized 

power in the presence of a third order nonlinearity. k 
is the nonlinear feedback parameter introduced by Sol- 
jacic et. aljl^]. n ^ l/KfF is an indicator of the extent 
of the field that is confined in the nonlinear region. 70 is 
the natural cavity linewidth. We obtained a Q = 80, 000 
and an on-resonance transmission T = 73% for the on- 
substrated and capped cavity. This corresponds to a 
Qsc = 500, 000, which is comparable with our previ- 
ously reported results for free-standing photonic crystal 
nanobeam cavity [isj. The demonstrated transmission is 
much higher than previous works with the same Q_25]. 

In summary, we proposed and demonstrated a deter- 
ministic design of the high Q PhC nanobeam cavities. 
Our cavities are ideally suited for the coupling to the 
feeding waveguides and allow for high transmission effi- 
ciency. This makes them ideal candidates for the real- 
ization of densely integrated photonic systems, and are 
suitable for applications ranging from optical intercon- 
nects to biochemical sensors. Finally, record high Q/V 
ratios of photonic crystal nanobeam cavities will enable 
further fundamental studies in spontaneous emission con- 
trol, nonlinear optics and quantum optics. 
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